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Abstract 

We present results from recent calculations of the QCD equation of state by the HotQCD Collab- 
oration and review the implications for hydrodynamic modeling. The equation of state of QCD 
at zero baryon density was calculated on a lattice of dimensions 32^ x 8 with mi - Q.I nis (corre- 
sponding to a pion mass of ~22Q MeV) using two improved staggered fermion actions, p4 and 
asqtad. Calculations were performed along lines of constant physics using more than lOOM cpu- 
hours on BG/L supercomputers at LLNL, NYBlue, and SDSC. We present parameterizations of 
the equation of state suitable for input into hydrodynamics models of heavy ion colUsions. 



1. Calculating the Equation of State with Lattice QCD 

Lattice Quantum Chromodynamics (LQCD) is the only viable method for calculating the 
equation of state for Quantum Chromodynamics in the vicinity of the transition between hadrons 
and the quark gluon plasma. As such, LQCD provides an essential input for modeling and 
thereby understanding the physics of relativistic heavy ion collisions. Using BG/L class super- 
computers at LLNL, BNL, and SDSC, the HotQCD collaboration has performed caculations of 
the equation of state and developed parameterizations suitable for input into hydrodynamic mod- 
els. The essential results for these calcualtions are summarized in this proceedings. Additional 
details and insights are found in Bazavov et al. 

As described in the seminal work of Wilson ||2| and Creutz |I3], the partition function for the 
QCD action can be calculated numerically on a discrete Euclidean lattice. In this form, the gluon 
action becomes a sum over all closed four-link loops (plaquettes) that can be evaluated with 
finite difference approximations. The full theory of QCD is recovered in the continuum limit, in 
which the lattice spacing is reduced to zero. Both the p4 and asqtad gauge actions are Symanzik 
improved, which means that they incorporate a six-link term that improves the convergence to 
0(a^). Both p4 and asqtad stagger the fermion spinors on odd-even lattice sites, which preserves 
a discrete chiral symmetry (i.e. exact chiral symmery is recovered in the continuum limit). Both 
actions also include O(a^) improvements, but they differ in their approach to further reducing 
cutoff effects. The p4 actions uses fat link smearing |4|, whereas the asqtad action uses tadpole 
coefficients [5|. 

The runs were performed along lines of constant physics, meaning that the bare quark masses 
are adjusted to maintain a constant value of the pion mass for a light quark mass set to one tenth 
the value of the strange quark mass. The temperature is determined by determining the distance 
ro at which the slope of the heavy quark potential equals l.65/rfy We use ro = 0.469(7) fm 
in this analysis Q. The thermodynamic quantities were calculated by averaging over approxi- 
mately 15k trajectories for 24 temperature values for p4, and 17 temperatures for asqtad, with a 
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Figure 1: Contributions the trace anomaly from gauge (left) and fermion (riglit). 



corresponding zero temperature runs of approximately 5k trajectories. A complete listing of run 
statistics, parameters, and thermodynamic quantities is given in the appendices of 1 1 1. 

The pressure and energy density are derived from the calculation of the interaction measure or 
trace anomaly, @^'^'{T)/T* - {e-3p)/T*, which measures the deviation from the conformal limit. 
The pressure is obtained by integrating above a temperature. To, for which the trace anomlay is 
close to zero. 
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The expression for the trace anomaly from the fermion and gluon components is, 

= -/;^7;,„<(2'?j/A(«A'A>/ + '?>.A<0'^),) , (2) 
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Here we use the notation A {X} - {X)(, - {X}^ to denote subtraction of the corresponding zero 
temperature values. /?„ is related to the derivative of the tadpole coefficient with respect to the 
gauge coupUng and is used in the asqtad calculations; this parameter is set to zero for p4. 



2. Equation of State Results 

The contributions from the gauge and fermion fields to the trace anomaly are shown sepa- 
rately in Figure [T] The dominant contribution is from the gluon action (left), although the cutoff 
effects appear to be larger for the fermion action (right). Note that this figure also includes results 
for Nr-^fi lattices. The full trace anomaly calculations for p4 and asqtad are shown in Figure[2] 
The two actions show good agreement in the low and high temperature regions, although the 
peak is lower for asqtad than for p4. Figure [2] also shows the parameterized fits of the lattice 
results, and fits to lattice results joined to hadron resonance gas calculations below 130 MeV. 

2 



0.2 0.4 0.6 0.8 1 



1.2 



Figure 2: Paramerized fits to tlie trace anomaly for p4 and asqtad calculations on lattices with = 8. 

3. Parameterizations for Hydrodynamics 

To facilitate the use of the HotQCD equation of state in hydrodynamic calculations we pro- 
vide a simple parameterization for the trace anomaly in Eq. Q. 











. P4 - 
1 asqtad ^ 
HRG • 




■'it 

■iff 

•w 
//W 








T [MeV] 




100 200 


300 400 


500 



e-3p 

J.4 



= 1- 



1 



d2 d4 

[1 +e(7--c,)/o)]2|^^ + f4 



(4) 



The parameters for these fits are listed in Table[T]for both the lattice fits and fits to the lattice and 
HRG results below 130 MeV. Additional details of the fitting procedure and systematic variations 
are given in Appendix C of Q. The pressure, energy density, and speed of sound squared are 



Data 


d2 [GeV]^ 


d4 [GeV]4 


ci [GeV] 


C2 [GeV] 




p4 

HRG+p4 


0.24(2) 
0.24(2) 


0.0054(17) 
0.0054(17) 


0.2038(6) 
0.2073(6) 


0.0136(4) 
0.0172(3) 


26.7/19 


asq 

HRG+asqtad 


0.312(5) 
0.312(5) 


0.00 
0.00 


0.2024(6) 
0.2048(6) 


0.0162(4) 
0.0188(4) 


34.4/14 



Table 1: Parameter values for fits of Eq. |4j to trace anomaly data for p4 and asqtad, and for data combined with HRG 
calculations. 



shown in Figure |3] The pressure was derived from Eq. |4] by numerical integration, beginning 
at a temperature of 50 MeV, the approximate minimum for this functional form. The shaded 
region for the energy density indicates the error associated with the linear vs. spline interpolation 
method described in |[T|, and the shaded region near the pressure shows the contribution from a 
hadron resonance gas calculation for the region below 100 MeV. Figure [3] also shows the speed 
of sound squared for a hadronic gas with a first order phase transition (double-dotted) as well as 
the default lattice inspired equation of state used by the 2D-I-1 viscous hydrodynamic model of 
Luzum and Romatschke [fl. 
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Figure 3: Energy density and pressure (left) and square of the speed of sound for parameterizations of the trace anomaly 
given by Eq.|4] 

We have presented new results on the equation of state of QCD with a strange quark mass 
chosen close to its physical value and two degenerate light quarks with one tenth of the strange 
quark mass. A comparison of calculations performed with the p4 and asqtad staggered fermion 
discretization schemes shows that both actions lead to a consistent picture for the temperature 
dependence of bulk thermodynamic observables. We have provided a simple parameterization 
for the equation of state suitable for insertion into hydrodynamic models of heavy ion collisions. 
Additional discussion of deconfinement and chiral symmetry restoration aspects of QCD ther- 
modynamics for the calculations on lattices with A^^ = 8 described herein may be found in [Ij. 
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